Front propagation in A^2A, A^3A process in Id: velocity, diffusion and velocity 

correlations. 
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We study front propagation in the reaction diffusion process {A A 2A, A -5- 3A} on a one 
dimensional (lei) lattice with hard core interaction between the particles. Using the leading particle 
picture, velocity of the front in the system is computed using different approximate methods, which 
is in good agreement with the simulation results. It is observed that in certain ranges of parameters, 
the front velocity varies as a power law of et, which is well captured by our approximate schemes. 
We also observe that the front dynamics exhibits temporal velocity correlations and these must be 
taken care of in order to find the exact estimates for the front diffusion coefficient. This correlation 
changes sign depending upon the sign of e t — D, where D is the bare diffusion coefficient of A 
particles. For et = D, the leading particle and thus the front moves like an uncorrelated random 
walker, which is explained through an exact analysis. 



I. INTRODUCTION 

Front propagation is an important field of study in 
nonequilibrium systems. We often encounter these prop- 
agating fronts separating different phases in physics, 
chemistry and biology Q. Here, in this work, we study 
the dynamics of the front in the reaction-diffusion sys- 
tem A — ► 2A, A — ► 3A in one dimensional lattice. At the 
macroscopic level, the mean field theory yields the fol- 
lowing partial differential equation for the coarse grained 
concentration p(x,t), 



dp 
dt 



D 



+ 2ep(l -p) + e t p(l - pj 2 



(1) 



where, D is the diffusion coefficient of the particle and e 
and et are the rates of single (A — > 2 A) and twin (A — » 
3A) offspring production respectively. Equation Q re- 
duces to the well known Fisher equation jfj when ej = 0, 
which models the reaction diffusion equation A — > 2A. 
The microscopic lattice model for A — > 2 A has been stud- 
ied extensively 00011. The mean field allows travel- 
ling wave solution of the form p{x, t) = <p(x — vt), where 
the velocity of an initially sharp front between p = 1 
( stable ) and p = ( unstable ) state approaches an 
asymptotic velocity Vq — 2^J (2e + e t )D. 



II. MODEL, FRONT VELOCITY AND 
DIFFUSION COEFFICIENT 



We consider a Id lattice composed of sites i = 1, 2...L. 
We start with the step function like distribution where, 
the left half is filled with A particles while the right half 
is empty. Each site i can either be empty or occupied 
by maximum one particle i.e. hard core exclusion is en- 
forced. We update the system random sequentially where 
L microscopic moves correspond to one Monte Carlo step 
(MCS). During each updating we randomly select a site 
and the particle at the site can undergo one of the fol- 
lowing three microscopic moves. 



(1) The particle can jump to neighbouring empty site 
with rate D, 

(2) The particle can give birth of one particle at either 
of the empty neighbouring sites with rate e, 

(3) The particle can generate two new particles at both 
the neighbouring sites provided both are empty with rate 
£*■ 

These processes are shown in the Fig. QJ. 
As time evolves, these stochastic moves result in the 
stochastic movement of the front. As has been argued 
in [5|, the front may be identified with the rightmost A 
particle. In this paper, we are interested in the dynamics 
of front whose evolution may be described by the follow- 
ing master equation 0. 



dP(X, t) 
dt 

DQ (X + 



= (e + D)P(X-l,t)- 
l,t)-(e + D)P(X,t) 



e t Q (X-l,t) + 
-(e t + D)Q (X,t) (2) 



Here, P(X, t) is the probability distribution of finding the 
front particle at position X at a time t and Qq (X, t) is 
the joint probability that the front is at X at time t and 
the site just behind it is empty. In Eq. (0), the first term 
corresponds to the forward hopping of the front particle 
from the position X — 1 to reach X due to birth of a 
single particle or diffusion. The second term corresponds 
to twin production at left and right neighbouring sites 
of site X — 1 and which results in the front moving from 
X — 1 to A. The third term corresponds to the backward 
hopping of the front particle from position X + 1 due 
to diffusion, provided site X is empty. Last two terms 
account for the possible jumps of the front particle from 
position A, which leads to front moving either at A — 1 
or A + 1. The dynamical properties of the front that 
we want to study are its velocity V and the diffusion 
coefficient Df which are defined as: 



V 



dt 



< X(t) > 



(3) 



2 



D 



f 



Id , , , 
2dt < ™~ 



< X(t) >f > 



(4) 



Where, < X(t) >= ^IP(X,i). Now we use the 
Eq. J2J and normalization = 1 and taking 

Qo{X) = (1 — px-i)P(X), where px-i is the probability 
that site X — 1 is occupied. Denoting px-i — Pi, we ob- 
tain the following expression for the asymptotic velocity 
and diffusion coefficient of the front. 



V = e + e t - pi(e t - D) 



D f = -{e + e t + 2D - Pl {e t 



D)} 



(5) 



(0) 



In order to obtain the velocity and diffusion coefficient 
we need to know pi, which is the density of site just 
behind the front. In for e t = 0, it was shown that 
front velocity approaches asymptotically the mean field 
value V = Vo in the limit — — » oo, while V = e + D in the 
opposite limit — — ► 0. But, when we are in between these 
two extreme limits we need to know pi and we expect 
similar features when e t ^ 0. There is no method to find 
this value exactly. Here we present some approximate 
analytic estimates for pi and hence the front velocity. 
In subsection A, we use fixed site representation method, 
where a truncated master equation is written in the frame 
moving with the front, as discussed in 0- In subsection 
B, we apply two particle representation scheme proposed 
by Kerstein £| while in C a mixed scheme is proposed 
which yields better results than either A or B. 



A. Fixed site representation 

This method has been proposed in |5| for the reaction 
diffusion process A <-> 2A. Here, we write a truncated 
master equation in the frame moving with the front. The 
simplest set of states is:{o»,»»}, which corresponds to 
the evolution of occupancy at a site just behind the front 
particle(Z = 1). Here the rightmost • in each state cor- 
responds to the front particle. These two states make 
transitions between each other due to the microscopic 
processes in the system as shown in Fig. . Considering 
all such transitions the evolution of probabilities of these 
two states are given by: 



dP(c 



dt 



dP(»> 

~dT 



= (2D-Dp 2 )P{.,)-{2Dp 2 + e(2 + p 2 ) 

+e t (l+p 2 (l-p 3 ))}P(o»), 
= {2Dp 2 + e(2 + p 2 )+e t (l + p 2 (l-p 3 )) 

}P(o.) - (2D - Dp 2 )P{..). (7) 



Here, pi is the density at the ith site behind the front and 
we have neglected the spatial density correlation between 
consecutive pairs of sites beyond the second site behind 
the front. Now, using Eq. Q and normalization P(o») + 
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FIG. 1: Microscopic moves, rightmost • represents the front, 
(a) Diffusion of the front particle to its right site leading to 
transition from •• to am with rate D. (b) Creation of one 
particle to the left of the front leads to transition from o« 
to •• with rate e. (c) om changes to •• due to creation of 
twins at both neighbouring sites of the front with rate et- (d) 
om — » ••, if the front takes diffusive move to its left and the 
second site behind the front is occupied. This occurs with rate 
Dp2, where p2 is the probability of occupancy at the second 
site behind the front. 



P(**) — 1, we obtain the following expression for p\\ 

2Dp 2 + e(2 + P2 ) + e t (l + p 2 (l - p 3 )) 
Dp 2 + 2D + 2e + ep 2 + e t (l + p 2 (l - p 3 )) 



Pi = 



(8) 



From Eq. © , we note that in order to find pi we need to 
know p 2 and p 3 . As a crude approximation if we assume 
that p 2 = p 3 = p b = 1, where p b is the bulk density, we 
get the following value of pi. 



Pi 



2D + 3e + e t 
3D + 3e + e t 



(9) 



Now using this approximation for pi in Eq. JSJ), we find 
the estimate for the velocity which is in reasonable agree- 
ment with the simulation, as shown in the Fig. (J2J . The 
estimate for V can be improved by including more sites 
in the truncated representation. For example, for I = 2 
we study the evolution of following set of four states: 
{o o •, o • •, • o • • •} and as expected we get improved 
results as shown in Fig. J3J|. Here, we notice that for 
larger values of et the simulation results show nice agree- 
ment with that of analytic results. However, as et de- 
creases and approaches zero, we see gradual departure 
of the simulation data from the analytic one. In fact, pi 
differs from the bulk density significantly with decreasing 
value e t as shown in the Fig. (jSJ. That is, the approxima- 
tion pi fa 1 holds better for larger values of et and hence 
we get better agreement with the simulation results. The 
estimate for the velocity can be further improved if we 
include states with larger number of sites. In Fig. Q 
we notice two interesting points: firstly, for D = e t , the 
theoretical result matches strikingly with the simulation 
result, secondly, we observe a power law dependence of 
the velocity on e t . In fact, the first point, can be shown 
to be exact by noting that when D = e t , the front ve- 
locity from Eq.JSJ is V = e + e t , which is independent of 
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FIG. 2: Transition between two particle states with rightmost 
• representing front, (a) Diffusive move of the front particle 
to its left leading to transition 101— >11 with rate D, (b) When 
the second particle behind the front jumps to the left, pro- 
vided it is empty, state changes from 101— >1001 with rate 
D(l — p), (c) Birth of a single particle by the second parti- 
cle to its left with rate e leads to transition 1001-^101, (d) 
1001— >101 if the second particle gives birth of two particles, 
provided the the site left to it is empty, with rate et(l — p). 



Pi- 



This equation specifies the dependence of p on the pa- 
rameters implicitly through dependence of Pq on e, et, D. 
Here, A is a free parameter to be evaluated as follows. 
Following Kerstein, using the ansatz P k = Pq(1 — Po) k 
and Eq. l(TT|) in Eq.JTJJl, we get the following quartic 
equation in Pq. 

e t AP 4 + (DX - e t - 2e t A)P 3 + (e + D + 2e t + e t X 

-DX)P$ + eP Q - 2e - e t = (12) 

In order to find Pq we need to fix the value of A. For large 
D and et = 0, it is known that the front particle moves 
with its mean field velocity 0- If we assume that this also 
happens when et ^ 0, then equating the mean field front 
velocity Vq — 2-\f (2e + et)D with that obtained from Eq. 
(J5J i.e. V ~ DPo when D is very large compared to other 

parameters, we get Pq = 2y/ 2e + £ * . Using this value of 

P Q in Eq. O we find A = 3/4 in the limit D -> oo. 
We solve the quartic equation l|12fl to get the value of 
Pi = Pq and hence the front velocity, as shown in the 
Figs. (@J and (0 and marked as 2P. 



B. Two particle representation 

In the following, we try to find the analytic estimates 
for pi using Kerstein's two particles representation^. In 
this representation each state of the system is defined by 
two rightmost particles and thus we have an infinite set 
of states: {11, 101, 1001, 10001, 100001, ....}. Here, the 
rightmost '1' represents the front particle while the left- 
most '1' is the second particle behind it and '0' stands 
for empty site. Let us denote by Pk the probability of 
two particle state with k empty sites between the leading 
particle and next particle behind it. These states form a 
closed set under transition due to microscopic processes. 
We have illustrated few transitions in the Fig. J2J). Con- 
sidering all such transitions and denoting the probability 
of occupancy of site just behind the second particle by p, 
we write the following rate equations for P^. 

dP 

= (e + 2D)Pi +e t (l-p)Pi + (2e + e t )(l-P ) 
-(2D-Dp)P , 

^ = (2D-Dp)P k _ 1 + {e + 2D + e t (l-p)}P k+1 
at 

- (4D- Dp + 3e + 2e t -e t p)Pk, k > 1. (10) 

In order to solve Ea. H10|) we need to specify the depen- 
dence of p on the parameters(e, et, D). Following Kerstein 
0, we write p = aPo — 6Pq and enforcing the condition 
that p = 1 when Pq = 1, we write the following expres- 
sion for p. 

P=(l + X)Pq - AP 2 (11) 



C. Mixed representation 



Since we are dealing with a multiparticle interacting 
system it is always desirable to include as many particles 
as possible while studying the evolution of the system. 
The simplest extension to the two particle representation 
is to study the evolution of the following set of states: 
{o o • o», • • o • ...}, where the rightmost • in 

each state denotes the front particle. Since in this rep- 
resentation, each state is characterized by two or three 
particles and hence we name it as mixed representation 
(MR). The rightmost • in each state denotes the front 
particle. When viewed in the frame moving with the 
front, each state contains the location of second parti- 
cle and the occupancy of the site just behind the second 
particle. We denote these states as (k, 0) or (k, 1), repre- 
senting the states having k empty sites between the front 
and the second particle and the site just after the second 
particle is empty or occupied respectively. For example, 
by (0, 0) we mean the state o • • and (0, 1) for the state 
• • •. These states are making transitions among each 
other due to the microscopic processes and form a closed 
set. We have shown some of the transitions in the Fig. 
©■ 

Assuming p as the density of site, which is, next nearest 
neighbour to the second particle, we write the following 
rate equation for the evolution of probabilities P(fc,0) 
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FIG. 3: Transitions between mixed particle states with the 
rightmost • representing front, (a) Diffusive move of the par- 
ticle to the right empty site with rate D.This leads to tran- 
sition from the state (1,1) to (0,0). (b) Birth of a new par- 
ticle on the right neighbouring empty site with rate e, which 
changes the state (1,0) to (0,1). (c) Transition from (1,0) to 
(1,1) with rate Dp, when the third particle jumps to the right 
neighbouring empty site, (d) (1,0)— >(0,1) when the second 
particle behind the front in (1,0) realization produces twins 
at the neighbouring empty sites with rate et- 
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FIG. 4: Percentage relative error in V , i.e. — , — 
(V s and V a representing simulation and analytic results re- 
spectively), is plotted against e t while keeping D=0.25 and 
e = 0.025 fixed. Here stars, open squares, open circles 
and filled circles correspond to the analytic estimates using 
I = 2(two states), I = 3 (four states), two particle(2P) and 
mixed representation (MR) respectively. Simulation profile 
and analytic profile using MR are essentially coincident and 
hence we notice almost zero error. Inset: log-log plot for the 
velocity versus et, the straight line shows power law V w e?' 25 . 



and p and then solve the equation: 



and P(k, 1), k = 0, 1, ..oo. 
dP(0,l) 



P(0,0) + P(0,1) = P 



(14) 



dt 



dP(0, 0) 
dt 



dP(k,l) 
Jt 



dP(k,0) 
dt 



{Dp + ep + 2e + e t p(l - p)}P(0, 0) 

(D + 2e + e t )P(l, 1) + (2c + 2e t )P(l, 0) 
e t {P(2, 0) + P(2, 1) + P(3, 0) + P(3, 1) 
(2D- Dp)P(0,l), 



+ ...} 



£>(l-p)P(0,l) + (D + e)P(l,l) 

(2D + e)P(l, 0) + 2e{P(2, 1) + P(2, 0) 
P(3, 0) + P(3, 1) + ...} - {2D + 2e + Dp 
ep + e t p(l- p)}P(0,0), 

DP(k - 1,1) + DpP(k- 1,0) 

{Dp + e + ep + e t p(l-p)}P{k,0) 

[D + e)P{k + 1, 1) + (e + e t )P(k + 1, 0) 

(4L> + 3e- Dp + e t )P(k,l), 

[D + D(l - p))P(k -1,0) + D(l - p)P{k, 1) 

DP(k + 1,1) + 2DP(k + 1,0)- {AD + 4e 
Dp + ep+2e t + e tP (l- p)}P(k,0). (13) 



In order to find Po, we need to solve the above set of 
coupled equations. However, one can find the analytic 
estimate for Po by solving rate equations for P(0,0) and 
P(0, 1) and assuming P(l, 1) = pPi,P(l,0) = (l-p)Pi. 

1 oo 

Using V^P(fc, i) — Pfe and V^Pfc = 1, we find steady 



i=0 



fe=0 



Following Kerstein (J], if we use the ansatz Pi = Po(l — 
Po), we get the following equation. 



aP (l-Po) + /3P +7-0 



(15) 



Where, 



state expression for P(0, 0) and P(0, 1) in terms of Pi 



a = (2e + e t + Dp-etp){3D + 2e + ep + e t p(l-p)} 
+ {2D -Dp- e){2D + 2e + ep + e t p(l - p)}, 

= {2D + 2e + ep + Dp + e t p(l-p)}(2D-Dp) 
-{Dp + 2e + ep + e t p(l - p)}(D - Dp), 

7 = e t (l- P ){3D + 2e + ep + e t p{l- p)} 

+2e(l - P ){2P + 2e + ep + e t p(l - p)}. (16) 

The Eq. (|15|l is in terms of two unknowns p and Po and 
hence we need to know p in order to find Pq . We specify 
the dependence of p on Po similar to what we did in 
the case of two particle representation. Using the results 
obtained for Pq = p\ in Eq.JSJ we find the estimates for 
V as shown in the Figs. I@J and © marked as MR. We 
observe good agreement with the simulation results. 



III. FRONT DIFFUSION COEFFICIENT 

In Fig. J7J), we have shown the simulation results of 
the front diffusion coefficient and compared it with the 
results obtained by using the mean field and simulation 
value of p\ in the equation©. Here, we notice the follow- 
ing interesting features: (1) when D — e t , the analytical 
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FIG. 5: Percentage relative error in V versus D, keeping e = 
0.05 and et = 0. The top data(filled squares) corresponds to 
ref.jfj, where in the frame moving with the front ,the evolution 
of particles at three sites behind the front was studied and 
assuming the fourth site at the bulk density. The middle 
data (open circles) corresponds to Kerstein p| two particle 
self consistent representation while the bottom data (filled 
circles) is the result from the mixed representation. 
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FIG. 7: Comparison of the front diffusion coefficient obtained 
analytically (from Eq. @ and using pi corresponding to I = 
2) with the simulation results for different values of et while 
keeping D = 0.25, e=0.025 fixed. We note that when e t =D = 
0.25, the simulation result matches with the analytic one. 



1 

0.9 
0.8 

0.7 f-/ 
0.6 
0.5 
0.4 





. • : 

.0"' 



: » • ! * ! 



10 

site(i) 



£,=0 
e t =0.05 
£,=0.15 
£,=0.25 
£,=0.35 



15 



20 



0.01 
0.005 


-0.005 
-0.01 
-0.015 



£,=0.05 — s 
£,=0.15 
£,=0.25 

£,=0.30 b 

£,=0.40 — . 



-6 
-9 
-12 

I 

10 

t 



15 



20 



FIG. 6: Density profile behind the front for different values 
of et while keeping D — 0.25 and e = 0.025 fixed. Here, we 
notice that as et decreases density profile curve shifts away 
from the bulk density level. 



FIG. 8: Simulation results for velocity correlation of the front 
with time t, for different values of et while keeping D = 0.25, 
e=0.025 fixed. We notice that when e t = 0.25, this correlation 
is zero. Inset: log-normal plot for C(t) versus t with D — 
0.05, 0.10, 0.30, 0.40 from top to bottom. 



value DJ" 1 " 1 matches well with the simulation result Df m . 
(2) when e t > D, the Df m > Df la (3) when e t < D, 
jjstm ^ D<j na . The origin of the above discrepancy be- 
tween the simulation and analytical results can be traced 
to the master equation @, where we have neglected the 
temporal velocity correlations. The expression for the 
asymptotic front diffusion coefficient with temporal cor- 
relations in velocity is given as: 

oo 

D f =D + J2C(t) (17) 

t=i 

where, Do is the front diffusion coefficient by neglecting 
correlations as given by Eq. © and C(t) is the temporal 



velocity correlation defined through, 
C(t) =< v(t')v(t' +t)>-< v{t') >< v(t' + t) >, 

(18) 

where, v(t) is the displacement of the front at time t. 
In Fig. ©, we have plotted the temporal velocity cor- 
relation C(t) for different values of e t . For e t > D, we 
observe positive correlation while for e t < D, it is nega- 
tive and for e t = D, C(t) seems to vanish for all t. Thus, 
e t = D is a special case, where the front particle moves 
like an uncorrelated random walker. 

In the following, we explicitly show that for the spe- 
cial case et = D, two consecutive steps of the leading 



6 



particle are uncorrelated, i.e., C(l) = in the steady 
state. Since at most two sites behind the front can be 
affected in two consecutive steps, we consider 4 states 
corresponding to I = 2, namely, {00f,0ff,101,fff} with 
the rightmost 'V representing the front. In order to find 
C(l) =< i>(i)t;(i+l) > - < v{t) >< v(t+l) >, we write 

< v(t)v(t + 1) >= R ++ - - +R , where Rij 

denotes the 'flux' R^ 1 + R^ 1 + R}® 1 + i? 111 for taking 

l j l j <-j l j 

two consecutive steps as i — +/— and j = +/—■ Here, 
for example, RP?}_ is the flux of two consecutive negative 
steps starting from the state 001. The only way it can 
occur is if the front particle takes two diffusive moves 
to the left and thus, = D 2 Pooi, where Pool is the 
steady state weight of the configuration 001. Considering 
all such two successive moves in the each state, we write 
the following expression for R-\ , R i- and R , 

R++ = D 2 + 2eD + e t D + e 2 + {e t D + e t e){P oi+Pioi}, 
R+- = D 2 , 

R-+ = (D 2 + De + De t )P Q0 i + (D 2 +De)P wl7 

= D 2 P 001 . (19) 

Now using Eas. H19|) . we have 

<v(t)v{t + l)> = e 2 + De t + 2eD + {e t D + e t e 
-D 2 -De)}P wl + {e t e-De}P 001 , (20) 

and similarly, 

< v(t) >=< v(t + 1) >= e + e t - {Pon + Pm}(e t - D). 

(21) 

Using Eas. lj2Q|l and Q21|l we can find the value of < 
v(t)v(t + 1) > or < v(t) > by finding the probabilities 
of different states, which is harder to compute exactly. 
However, when e t = D, we find that < v(t)v(t + 1) > 
— < v(t) >< v(t + 1) > is independent of all the proba- 
bilities and is equal to zero, i. e., two successive steps are 
uncorrelated as observed in the simulation result (FigJSJ. 
We also note that for this special case the above analy- 



sis does not involve any approximation i.e. it is exact. 
We also notice that when e t ^ D, < v(t)v(t + 1) >/< 
v(t) >< v(t + 1) >, i.e., the front motion is correlated. 
Preliminary fits suggest that the the temporal velocity 
correlation has the form, C(t) ~ t a e~^ 1 . 



IV. CONCLUSION 

In the present paper, we studied the reaction diffusion 
system A — > 2A, A — > 3A in one dimension. Treating the 
rightmost occupied site as a front we compute the front 
velocity analytically using different approximate meth- 
ods. In fixed site representation one can systematically 
improve upon the estimate by studying the evolution of 
particles at larger number of sites behind the front. The 
results from two particle representation and mixed repre- 
sentation show excellent agreement with the simulation 
results. We also observed that the velocity depends on e t 
as a power law. As far as the computation of front dif- 
fusion coefficient is concerned, we notice that one needs 
to take into account the temporal velocity correlation. 
In fact the observed temporal correlations in the front 
dynamics changes sign with sign of e t — D. For e t > D 
or et < D, front moves like a positively or negatively 
correlated random walk and for e t = £>, the temporal 
correlations in steps vanish and the front particle moves 
like a simple uncorrelated random walker. An interest- 
ing generalization of the process would be to include the 
annihilation of particles as well i.e., 2 A — > A with rate 
W. The case e t = and W — D is exactly solvable [|| 
and both temporal and spatial correlations vanish. For 
non-zero D, e t and W, simulations show that temporal 
correlations vanish on the plane D = e t + W, although 
spatial correlations do not. Computation of the previous 
section can be extended to this case to show that consec- 
utive steps are temporally uncorrelated. A general proof 
to show that C (t) vanishes for all t remains an interesting 
open problem. 
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